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We formulated an effective theory for a single interlayer exciton in a bilayer quantum antiferro- 
magnet, in the limit that the ho Ion and doublon are strongly bound onto one interlayer rung by the 
Coulomb force. Upon using a rung linear spin wave approximation of the bilayer Heisenberg model, 
we calculated the spectral function of the exciton for a wide range of the interlayer Heisenberg 
coupling a = J \/ 'J r z. In the disordered phase at large a, a coherent quasiparticle peak appears 
representing free motion of the exciton in a spin singlet background. In the Neel phase, which 
applies to more realistic model parameters, a ladder spectrum arises due to Ising confinement of the 
exciton. The exciton spectrum is visible in measurements of the dielectric function, such as c-axis 
optical conductivity measurements. 
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I. INTRODUCTION 



An exciton is the bound state of an electron and a hole, 
and considering their bosonic character the question im- 
mediately arises whether they can condense into an exci- 
ton Bose condensate^. The quest for such exciton super- 
fluidity has, over the past decade, increasingly focussed 
its attention to layered structures where one layer con- 
tains holes and the other layer contains electrons 2 . The 
Coulomb attraction between the electrons and holes then 
allows for the formation of so-called interlayer excitons. 
In 2004 a condensate of interlayer excitons was success- 
fully created in a heterostructure of two 2DEGs under 
the application of a perpendicular magnetic fielcP. Since 
then many other candidate materials were suggested that 
should support interlayer exciton condensation in the ab- 
sence of magnetic fields, such as graphene^"^ or topolog- 
ical insulators^. One class of candidate materials has 
not been considered yet, namely the Mott insulators 9 . 
The strong interactions between electrons make these 
materials currently one of the most fascinating and the 
least understood solid state compounds. When making 
heterostructures of p and n-doped quasi-two-dimensional 
Cu02 layers one expects the formation of interlayer exci- 
tons, and these excitons will interact strongly with mag- 
netic excitations, possibly leading to unexpected dynam- 
ics. To explore all these unexpected dynamics of the exci- 
tons in the strongly correlated system such as the exciton 
condensation, understanding the dynamics of single ex- 
citon will be the first step. 

Heterostructures of p and n-doped cuprates can be typ- 
ically described by a strongly correlated model: the bi- 
layer t — J model, which is extended from two single-band 
t — J model for each layer with coupling terms between 
the layers as following: 



H bt _j = H t + Hj + H v (1) 



where H t is the hopping of electrons in each layer 

H t = -t e 4lo C Jl* + h - C - ( 2 ) 

(ij)o;l 

and Hj is the bilayer Heisenberg model describing the 
undoped Mott insulating state 

Hj = J ^ S U ' S jl + J ± ^ S H ' S i2- (3) 

Here cu a and Su denotes the electron and spin operators 
respectively on site i in layer I = 1,2. The Heisenberg 
Hj is antiferromagnetic with J > and J±_ > 0. The 
last term Hy in ([TJ is the Coulomb attraction between 
a vacant site (holon) and double-occupied site (doublon) 
in the same rung, described by 

H v = ^2nnni2 (4) 

i 

which is the force required to form an exciton in the same 
rung. Without loss of generality, we assume that layer '1' 
has n-type doping with the constraint ^ a c\ la cn a > 1 
and layer '2' has p-type with the constraint ^ a c| lf7 Qi a < 
1. 

In this paper we will present a theoretical framework 
describing the dynamical properties of a single exciton in 
a strongly correlated bilayer described by ([!]), following 
our previous shorter publication on this topic 10 . The 
binding of the holon and the doublon is determined by 
the interlayer Coulomb repulsion and we will focus on 
the strong coupling limit (V > t). This implies that the 
exciton is formed by the holon and doublon on the same 
rung, as is shown in Figure [l] 

Understanding of the bilayer Heisenberg model will be 
an important step towards analysing the dynamics of 
a single exciton. The ground state and excitations of 
the bilayer Heisenberg Hamiltonian have been studied 
quite extensively using Quantum Monte Carlo (QMC) 
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FIG. 1: Naive real space picture of an exciton in a strongly 
correlated material, as viewed from the side. Two square 
lattices (blue balls) are placed on top of each other. The 
red arrows denote the spin ordering, which forms a perfect 
Neel state. The exciton consists of a bound pair of a double 
occupied and a vacant site on an interlayer rung. The energy 
required to break this doublon-holon pair is V. The magnetic 
ordering is governed by the in-plane Heisenberg J and the 
interlayer Jj_, as described by the Hamiltonian 
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FIG. 2: Zero temperature phase diagram of the bilayer 
Heisenberg model as a function of interlayer coupling strength 
a = on the horizontal axis. At a critical value a c a quan- 
tum phase transition exists from the antiferromagnetic to the 
singlet phase. The vertical axis shows the Neel order parame- 
ter signaling antiferromagnetism. Note that even at a — the 
Neel order parameter is reduced from the mean field value \ 
to approximately 0.3 due to spin flip interactions. (Adapted 
from RefP.) 



method^ 11 * 12 *, dimer expansiomP^^ and the closely re- 
lated bond operator theor^^l, the nonlinear sigma 
model^U^ anc [ S pi n wa ve theorjPDHSD All results indi- 
cate a 0(3) universality class quantum phase transition 
at a critical value of J_\_/ J from an antiferromagnetically 
ordered to a disordered state, see Figure [2j A naive mean 
field picture of the antiferromagnetic ground state is pro- 
vided by the Neel state, in which each of the sublattices 
are occupied by either spin up or spin down electrons as 
shown in Figure [I] However, the exact ground state is 
scrambled up by spin flip interactions reducing the Neel 
order parameter to about 60% of its mean field valued 
A finite interlayer coupling J± influences the antiferro- 
magnetic order. In the limit of infinite Jj_, the electrons 
on each interlayer rung tend to form singlets destroying 
the antiferromagnetic order. 

Standard spin wave theories however cannot account 
for the critical value of J±/J ~ 2.5 found in QMC and 
series expansion studies. This discrepancy between nu- 
merical results and the spin wave theory has a physical 



FIG. 3: Exciton motion in a naive real space picture. In 
a perfect Neel state, the motion of an exciton (with respect 
to the situation in Figure [T]) causes a mismatch in the spin 
ordering. The kinetic energy gained by moving the exciton is 
proportional to the energies of the doublon t e and holon tu 
divided by the exciton binding energy V. 



origin. Chubukov and Mori 25 pointed out that standard 
spin wave theories do not take into account the longi- 
tudinal (that is, the interlayer) spin modes. By tak- 
ing into account those longitudinal spin waves one can 
derive analytically the right phase diagram^. Another 
correct method is to introduce an auxiliary interaction 
which takes care of the hard-core constraint on the spin 
modest 

If one wants to study the doped bilayer antiferromag- 
net however, one needs explicit expressions of how a mov- 
ing dopant (be it a hole, electron or exciton) interacts 
with the spin excitations. Even though the Neel state is 
just an approximation to the antiferromagnetic ground 
state, it provides an intuitive explanation of the major 
role spins play in the dynamics of any dopant. As can 
be seen in Figure [3J a moving exciton causes a mismatch 
in the previously perfect Neel state. Consequently, the 
motion of an exciton is greatly hindered and a full under- 
standing of possible spin wave interactions is needed to 
describe the exciton dynamics. This is of course similar 
to the motion of a single hole in a single Mott insula- 
tor layeiJ^HMl. Vojta and BeckeJ^ have computed the 
spectral function of a single hole in the Heisenberg bi- 
layer. A rung linear spin wave approximation 26 is needed 
to obtain the expressions for the spin waves in terms of 
single site spin operators. Summarizing, we will formu- 
late first an effective exciton t — J model from the bi- 
layer t — J model in the limit of strong Coulomb at- 
traction in section |TTJ In order to find the interaction 
coefficients between excitons and spin excitations we will 
construct a spin wave theory of the bilayer Heisenberg 



model in section III Based on these two developments, 
we can compute the exciton spectral function using the 
self-consistent Born approximation in section [TV) Finally, 
we connect the exciton spectral function to measurable 
quantities in section [V] 



II. THE BILAYER EXCITON t - J MODEL 

The bilayer t — J model describes generally the 
p/n-doped bilayer antiferromagnet. The behavior of a 
bound exciton however depends on the magnitude of the 
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FIG. 4: In weak coupling the spectrum of an exciton is ob- 
tained by the ladder diagram approximation from the spec- 
trum of the single doped hole. The Xo an d Xo are respectively 
the imaginary and real part of the bare exciton susceptibility. 
The x and x are the imaginary and real part of the full 
exciton susceptibility obtained in the ladder diagram approx- 
imation |5j). Besides the continuous particle-hole spectrum 
above the gap, there can only be a single exciton peak deter- 
mined by Vxo = 1 in the weak coupling limit. 



Coulomb force V in Hy, equation Q. If the Coulomb 
repulsion is relatively weak, the motion of holons and 
doublons will be relatively independent with each other 
and the Hy can be treated as a perturbation on H t + 
Hj. The full exciton-susceptibility xi 00 ) can be obtained 
from the bare susceptibility Xo(^) m the absence of the 
Coulomb force using the ladder diagram approximation, 



(5) 



Since the undoped state is a Mott insulator, there is a 
gap in the imaginary part of the bare susceptibility Xo- 
Above this gap there is an onset of the particle-hole con- 
tinuum. In the ladder diagram approximation, there can 
only be a single delta function peak in the full suscepti- 
bility at Vx'o — 1 signaling the formation of an exciton. 
We conclude that in the weak coupling limit no special 
exciton features other than a single delta function peak 
can appear in the gap. Following our expectation that 
realistic materials are in fact in the strong coupling limit, 
as explained in section [VJ we will henceforth focus our 
attention to the strong coupling limit. 

In the strongly coupling limit (V ^> t), the hopping 
term H t can be treated as a perturbation on the unper- 
turbed Hy using the perturbation method developed by 
KatcW In this method, one considers first an exact solv- 
able part of the Hamiltonian, in this case the interlayer 
Coulomb interaction Hy. It has the eigenvalues 



V(N - N + N) = Eq + VN 



(6) 



where N is the total number of sites, Nq is the number of 
dopants per layer and N is the number of double occupied 
sites that do not lie above a vacant site. It is clear that 
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FIG. 5: The motion of the composite exciton can be re- 
lated to the motion of its constituents via Kato's perturbation 
method. In this method a virtual intermediate breakup of the 
exciton is in between the initial state (Figure [TJ and the final 
state (Figure The kinetic energy of the exciton is therefore 
the product of the kinetic energies of the holon and doublon 
divided by the energy of this virtual state, t ex = t 2 /V . 



the ground state of Hy is given by the state where all 
double occupied and vacant sites lie above each other, 
as depicted in Figure [I] As mentioned before an exciton 
consists of a double occupied and a vacant site bound on 
top of each other. Consequently, the ground state of Hy 
is the state where all dopants are bound into excitons. 

The essence of Kato's perturbation method is that we 
now forbid all states with higher Hy eigenvalues. In our 
model, this implies that we forbid states such as the one 
depicted in Figure [5] where the double occupied site is 
not on top of the vacant site. In zeroeth order, hopping 
of electrons is forbidden since that would break up an 
exciton state. Therefore the zeroeth order Hamiltonian 
only contains Heisenberg terms H^ = Hj. 

In second order processes are allowed that virtually 
break up excitons, but end up with only bound excitons. 
The corresponding effective Hamiltonian is given by 



2V 



P e (H t )(l-P e )(H t )P e 



(7) 



where P e is the operator that projects out states with 
unbound dopants. As can be verified from Figure [5] this 
process allows the hopping of excitons by virtually break- 
ing the dopants apart. If we define the exciton operator 
in terms of electron creation operators 



E 



t -J 



,t C L|(^ Pip) 



the exciton hopping process can be formulated as 



H 



t,ex 



teth 



E 4 

<ij>aa' 



Ana' C lpa C 3P(r C jn(T' 



(8) 



(9) 



Note that in this Hamiltonian, no break-up of the exciton 
is required. The virtual process as described before only 
enabled us to relate the single layer kinetic energies to 
the bilayer exciton kinetic energy, 



t 



t e th 
V ' 



(10) 



Here t e is the hopping energy for a single electron, th 
the hopping energy for a single hole and t is the hopping 



4 



energy for a bound exciton. In addition to this hopping 
process there are also second order processes that equal 
a shift in chemical potential of the excitons. In the limit 
that we are interested in, that of a single exciton, we 
neglect chemical potential terms. 

In conclusion, we formulated a model for the strong 
coupling limit of Hy that describes the motion of bound 
excitons in a Mott insulator double layer. The corre- 
sponding Hamiltonian is 



d — 2 , 



H = H, 



Hj 



(11) 



We will refer to this model as the exciton t — J model 



A. The singlet-triplet basis 

The hopping term ([9| represents an exciton Ei on site 
i swapping places with the spin background Cj pa Cj na > on 
site j. This Hamiltonian is in the electron Fock state rep- 
resentation with the background determined by the bi- 
layer Heisenberg model (|3| . Historically the spin singlet- 
triplet basis turned out to be convenient in treating the 
bilayer Heisenberg model, and consequently we will apply 
this representation also to the hopping term 

Unlike the fermionic holes in the single layer case, the 
exciton is composed of a fermionic doublon and holon in 
the same rung, and hence is a bosonic particle. The local 
Hilbert space on each interlayer rung is five dimensional 
with basis in terms of five hard-core bosons as one inter- 
layer exciton state \E)i and four different spin states. In 
the single-triplet basis, which is valid for both the doped 
and undoped case, we can introduce the four hard core- 
boson as one singlet state and three triplet states: 

loo)* = ^(4A-4 u 4 t )i°> (12) 

|io>i = ^(4it c m + 4u4t)l°) (1 3 ) 

ii i>< = 4 t 4 t |o> (14) 

|1 -1); = 4 t 4 t |0>. (15) 
Then the hopping term [9] can be reexpressed as: 

H t = -t J2 \ E i) (lOOMOOl. +^llmMlmU (^|. 

<ij> \ m J 

(16) 
(17) 

(18) 



We can introduce the the total spin operator 

Si = S^i + Si2 

and the spin difference operator 

S = Sn — Se- 



this reads 

si = |1 1)(1 1 

s? 
s+ 



1-1)(1-1| (19) 

>/2(|ll>(10| + |10)<l -1|) (20) 

-|0 0)(10|-|10)(0 0| (21) 

>/2(|ll>(0 0|-|0 0)<l -1|). (22) 



In general, we see that the operator Si conserves the to- 
tal onsite spin, while S always changes the total spin 
number s by a unit. The ^-components of the spin oper- 
ators do not change the magnetic number m, while the 
^-components of the spin operators change the magnetic 
number by a unit. The bilayer Heisenberg model is now 
written as 



H 



•s,- 



Si ■ Sj 



<ij> 



4 



£ 0? - ) • 



(23) 

In conclusion, we formulated the exciton t-J model in the 
singlet-triplet basis which will be a starting point to solve 
the dynamics of the single exciton. 



B. Sign problem 

Notice also that the Hilbert space no longer contains 
fermionic degrees of freedom. The question is whether 
the disappearance of the fermionic structure also leads 
to the disappearance of the fermionic sign structure, 
which causes so much difficulties in the single layer t — J 
modeP 5 -! 

The sign structure can be investigated as follows. Re- 
member that at half-filling the fermionic signs in the stan- 
dard t — J model on a bipartite lattice can be removed 
by a Marshall sign transformation^. Upon doping, signs 
reappear whenever a hole is exchanged with (for exam- 
ple) a down spin. Which matrix elements of the Hamil- 
tonian become positive (and thus create a minus sign in 
the path integral loop expansion) depends on the specific 
basis and on the specific Marshall sign transformation. 

For the double layer exciton model, define a spin basis 
state with a built-in Marshall sign transformation of the 
form (compare to Ref.^) 



= (-i) 



Itlt 
I I 



(24) 



where N^ n is the number of down spins on the A sub- 
lattice in the n-layer and similary we define N^ p . With 
these basis states the Heisenberg terms are sign-free and 
the only positive matrix elements come from the ex- 
change of an exciton with a m = ±1 triplet. 

We conclude that, even though the model is purely 
bosonic, the exciton t — J model is not sign-free and it is 
not possible to remove this sign structure using a Mar- 
shall or similar transformation. 52 However, as will be fur- 



Explicitly in terms of singlet and triplet rung states for 



ther elaborated upon in section IV , for both the antiferro- 
magnetic and singlet ground states these signs do cancel 
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out. Therefore for such ordered bilayers the problem of 
exciton motion turns out to be effectively bosonic. 



III. UNDOPED CASE: THE BILAYER 
HEISENBERG MODEL 

Before considering the dynamics of the exciton and ex- 
pressing the interaction between the exciton and the spin 
background, we need to derive a spin wave theory for 
the bilayer Heisenberg model. Similar to the traditional 
Holstein-Primakoff spin-wave theory, we need a classical 
reference state, i.e. the mean field ground state of the bi- 
layer Heisenberg model, then develop the linear order for 
the spin wave theory from the mean field ground state. 
The method we present here is similar to the one pre- 
sented in^l 



In the limit of large S these excitations tend to become 
purely bosonic. We will take the mean field ground state 
of (25) and these three excitations as the starting point 



for the linear spin wave theory. 

Finally, we must mention the obvious flaw in the above 
reasoning. Where we criticized earlier spin wave theo- 
ries because they predicted the wrong critical value of 
J±/Jz, we now apparently adopt such a 'wrong' theory 
since (26) predicts a c = 1 for S = \\ Nevertheless, as we 
show in appendix |C| concerning S = |, the presence of 
spin waves changestne ground state energy which makes 
the disordered state more favorable even below the mean 
field critical (jj) c calculated above. Hence, due to cor- 
rectly taken the ground state energy shifts into account, 
one finds the accurate critical value for a consistent with 
numerical calculations. 



A. Mean field ground state 



The singlet-triplet basis ( 23 ) of the bilayer Heisenberg 
model is convenient for mean field theory. Mean field the- 
ory tells us that for large ratio J±/J the ground state is 
the singlet configuration |0 0). For small J_l/J, we expect 
antiferromagnetic ordering, which amounts to staggered 
condensation of S z . By setting (S z ) = {—l) l m we obtain 
a mean field Hamiltonian 

h mf = £ [1 Jz ~2 + J± ( 5f _ ^ _ 1 j^^ys- 

(25) 

which has a ordered-disordered transition point at 



70 =5 S < S + 1 » 



(26) 



where S is the magnitude of spin of the spin operator on 
each site. A proof of this result can be found in appendix 

The basic idea of a spin wave theor}E^^ is to start 
from this semiclassical (mean field) ground state and de- 
scribe the local excitations with respect to this ground 
state. One can immediately infer why the Holstein- 
Primakoff or Schwinger approach to spin wave theories 
fails for the bilayer Heisenberg model. First, the mean 
field ground state is no longer a Neel state for finite 
a. Secondly, where Holstein-Primakoff describes one and 
Schwinger describes two onsite spin excitations, the bi- 
layer Heisenberg has in fact three types of excitations. 
This has been pointed out by Chubukov and Morr^, who 
called the 'third' excitation the longitudinal mode. 

Here we want to point out that due to the local Hilbert 
space and the mean field ground state as described by 
(25) we can 'reach' all states in the local Hilbert space 
with three types of excitations: a longitudinal which 
keeps the magnetic number m constant, and transversal 
b± who change the magnetic number m by either ±1. 



B. Spin wave theory 

We will now construct explicitly the spin wave theory 
described above for S = \. First, one needs to find the 
ground state following equation (25). In the S = \ case, 



this amounts to a competition between the singlet state 
\s — 0, m — 0) and the triplet \s = 1, m = 0). The 
mean field ground state on each rung is given by a linear 
superposition of those two, 



\G)i = rji cosx|0 0)i - sinx|l 0)i 



(27) 



which interpolates between the Neel state (x = tt/4) and 
the singlet state (x = 0). The onset of antiferromag- 
netic order can thus be viewed as the condensation of the 
triplet state in a singlet background.^ With rji = ( — l) z 
alternating we have introduced a sign change between the 
two sublattices A and B. The angle x will t> e determined 
later by self-consistency conditions. 

The three operators that describe excitations with re- 
spect to the ground state are 



e\ = (rji sin x\0 0)i 



•cos X |10)0(G|, 



bl = |1 -l)i(GU. 



(28) 
(29) 
(30) 



The e-operators will later turn out to represent the lon- 
gitudinal spin waves, whereas the 6-operators represent 
the two possible transversal spin waves. 

The bilayer Heisenberg model can be rewritten in 
terms of these operators. For completeness we include 
the parameter A that enables a comparison with the Ising 
limit (A = 0) with the Heisenberg limit (A = 1), 

Si • S 2 = SIS Z 2 + ^A(S+S 2 - + S^S+). (31) 

Given this, we can explicitly write down the spin opera- 
tors in terms of the new e and b operators, as is done in 
appendix [Bl 
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From the requirement that the Hamiltonian does not 
contain terms linear in spin wave operators we obtain the 
self-consistent mean field condition for the ground state 
angle x, 



(cos 2\ — ct\) sin 2\ — 



(32) 



which has two possible solutions. Either \ — 0? which 
corresponds to a singlet ground state configuration, the 
disordered phase. If cos 2\ = a A, there exists an anti- 
ferromagnetic ordered phase. These are indeed the two 
phases represented in Figure |2j Which of the two solu- 
tions ought to be chosen, depends on the ground state en- 
ergy competition. In appendix [C] we compare the ground 
state energy of both phases, from which we can deduce 
that the critical p oint li es at a c « 0.6, consistent with 
numerical literatur e 1 1 * 12 l 

The dispersion of the spin wave excitations can be 
found when if we consider only the quadratic terms in 
the Hamiltonian. This is called the 'linear' spin wave ap- 
proximation, and it amounts to neglecting the cubic and 
quart ic interaction terms. First take a Fourier transform 
of the spin wave operators 



ot - 



2 ^ +- 

N Z^ e ko 

k 



ik-Tn 



(33) 



where the sum over k runs over the 2/N momentum 
points in the domain [— 7r, tt] x [— 7r, tt] and cr = A,B 
represents the sublattice index. A similar definition is 
used for the 6-operators. 

Upon Fourier transformation, we can decouple the spin 
waves from the two sublattices A and B by introducing 



J 

^k,p 



1 

7i 



(4a+p4b) 



(34) 



where p = ± stand for the phase of the spin mode. Modes 
with p = — 1 are out-of-phase and have the same disper- 
sion as the in-phase p = 1 modes but shifted over the 
antiferromagnetic wavevector Q = (7r,7r). Again similar 
considerations hold for the b operators. 

Next we perform the Bogolyubov transformation on 
the magnetic excitations, 



ot - 



cosh<pfc,pCfe, p + sinh<pfc, P C-fc,; 



(35) 



b{ iPi+ = cosh0 k , p al p + smh6 k , P f3-k,p 



(36) 

8wh.0 k ,pOt-k, p (37) 



The corresponding transformation angles are set by the 
requirement that the Hamiltonian becomes diagonal in 
the new operators £ (the longitudinal spin wave) and 
a,/3 (the transversal spin wave). In doing so, we intro- 
duced the 'ideal' spin wave approximation in which we 
assume that the spin wave operators obey bosonic com- 
mutation relations 40 . This assumption is exact in the 
large S limit. For S = \ this approximation turns out 
to work extremely welP^, since the corrections to the 



bosonic commutation relations are expressed as higher 
order spin-wave interactions. The Bogolyubov angles are 
therefore given by 



tanh 2(fk,p 



-p\ cos 2 



tanh 20 



k,p 



sin 2 2% + Aa cos 2\ — p\ cos 2 2x7& 
P^lk 



, (38) 



sin 2\ + (1 + A) a cos 2 x — pA cos 2%7i 



The factor j k encodes for the lattice structure, and it 
equals for a square lattice 



Ik 



- ^2 e%k ' 5 = \ ( cos k x + cos ky) 



(40) 



where the sum runs over all nearest neighbor lattice sites 
5. The Bogolyuobov angles still depend on x> which 
characterizes the ground state. In the antiferromagnetic 
phase cos2x = Aa and for the Heisenberg limit A = 1 
these angles reduce to 



tanh 2(fik,p 
tanh 20k, p 



-pct 2 lk 

- va 2 lk ' 
Vlk 



1 + a - paj k 



(41) 
(42) 



We can distinguish between the longitudinal and 
transversal spin excitations, with their dispersions given 

by 



"k,p 



Jzy/T 



-Jz 



v /(l + a(l-p 7 fc)) 2 -7| 



(43) 
(44) 



The longitudinal spin wave is gapped and in the limit 
where the layers are decoupled (a = 0) completely non- 
dispersive, while the transversal spin wave is always lin- 
ear for small momentum k. This type of spectrum is 
similar to a phonon spectrum, which contains a linear k- 
dependent acoustic mode and a gapped flat optical mode. 
This correspondence between spin waves and phonons 
enables us to use techniques from electron-phonon inter- 
action studies for the exciton-spin wave interactions. 

On the other hand, in the singlet phase (a > 1) one 
has trivially three identical triplet spin excitations. The 
Bogolyubov angles are given by 



tanh2(/?^p = — tanh 20 k , p 



-Plk 



2a - pj k 

and the dispersion of the triplet spin waves is 
£k, P = Jz^/a(a- pj k ). 



(45) 



(46) 



These dispersions cor respond to earlier numerical and se- 
ries expansions results 13 * 14 * 25 * 27 ^ . In fact, these results are 
exactly equal to the dispersions obtained from the non- 
linear sigma modePH 

The above derivation adds to earlier studies of the bi- 
layer Heisenberg model in that we now found explicit ex- 
pressions of how the spin waves are related to local spin 
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FIG. 6: Dispersion of the bilayer Heisenberg spin waves for 
different values of a. The top row has a = 0.04 and a = 0.4, 
the bottom row a = 0.9 and a — 1.1. In the antiferromagnetic 
phase (first three pictures) there is a clear distinction between 
the longitudinal spin waves (in green) and the transversal spin 
waves (in blue and red). The first is gapped, whilst the latter 
is zero at either k = (0,0) or (7r,7r) with a linear energy- 
momentum dependence. In the singlet phase, all spin waves 
are gapped triplet excitations. 



flips, equations (35)- (39). This microscopic understand- 
ing of the magnetic excitations of the system enables us 
in the next section to derive exactly how magnetic inter- 
actions influence the dynamics of excitons. 



IV. A SINGLE EXCITON IN A CORRELATED 
BILAYER 

As was pointed out in section [TTJ the exciton t — J 
model is still troubled by the sign problem even though 
it is purely bosonic. The sign-problem makes it difficult 
to say anything conclusive for systems with a finite den- 
sity of excitons. Doping the single layer t — J model 
leads to similar loss of theoretical control, and is the 
consequence of the fact that the magnetic ground state 
changes rapidly with doping. However, we can derive the 
dynamics of a single exciton in the undoped bilayer. In 
the thermodynamic limit a single exciton will not change 
the ground state. Following the exciton hopping Hamil- 
tonian ([9| we can express the dynamics of the exciton 
upon interaction with the spin wave modes. A single 
exciton can be physically realized by either exciting a in- 
terlayer charge-transfer exciton in the undoped bilayer, 
or by infinitesimal small chemical doping of layered struc- 
tures. 



From a theoretical side, the spin wave we derived in the 
last section HH can be used to constructed the effective 
theory of the single exciton and apply the self-consistent 
Born approximation. Similar to the single layer case^, 
we consider the mean field state \G) as the vacuum state 
and it is straightforward to derive the effective theory for 
single exciton as: 



H 



+ sin 2x(eJ 



cos2x(l - e\ej) 



(47) 



The dynamics of a single exciton are contained in the 
dressed Greens function, formally written as 



G*(k,u>) = (il>o\E kiP 



H 



HJ^o) (48) 



where E s k is the Fourier transformed exciton creation 
operator, and p indicates the same phase index as used 
for the spin waves in equation (34). The |^o) denotes 



the ground state after the spin wave theory has been 
applied. The Greens function cannot be solved exactly 
and one needs to write out a diagrammatic expansion in 
the parameter t/J. For this purpose, we have derived the 
corresponding Feynman rules of the exciton t — J model 
in appendix [P] 

Using Dyson's equation one can rephrase the diagram- 
matic expansion in terms of the self-energy £ p (/c, uj) such 
that 



G p (k,u) 



uj - eg(fc) - YP(k,u) + : 



(49) 



where e p (k) is the dispersion in the absence of spin ex- 
citations for the exciton with phase p. The self-energy 
can be computed by summing all one-particle irreducible 
Feynman diagrams. The degree to which exciton motion 
contains a free part grows with a, and indeed the free 
dispersion is 



eg(fc) = p zt cos 2x 7fc 



(50) 



where cos 2% equals aX in the antiferromagnetic phase 
and equals 1 in the singlet phase. 

As we noted before, the spin wave spectrum resembles 
a phonon spectrum. Hence we can compute the exciton 
self-energy using the Self- Consistent Born Approxima- 
tion (SCBA) 28 29 , an approximation scheme developed 
for electron-phonon interactions but subsequently suc- 
cessfully applied to the single layer t — J model. 

The SCBA is based on two assumptions: 1) that one 
can neglect vertex corrections and 2) one uses only the 
bare spin wave propagators. The first assumption is mo- 
tivated by an extension of Migdal's theorem^, the second 
by the linear spin wave approximation. Consequently, all 
remaining diagrams are of the 'rainbow' type which can 
be summed over using a self-consistent equation. For the 
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FIG. 7: Feynman diagram representation of the Self- 
Consistent Born Approximation (SCBA). The self-energy of 
the exciton depends self-consistently on 'rainbow' diagrams 
where it emits and absorbs either one or two spin waves. The 
left two diagrams contain interaction with the longitudinal 
spin wave, the right diagram contains the interaction with 
the transversal spin waves. The coloring is the same as in 
the dispersion diagram, Figure [6] Vertex corrections are ne- 
glected in the SCBA. 



exciton t — J model, the SCBA amounts to computing 
the self-energy for the in-phase exciton, as shown dia- 
grammatically in Figure [7] Analytically we write for the 
in-phase exciton self-energy, 



E+(Jfe,o;) 



zH 2 
N 



zH 2 



sin 2 2x ^2 (^ fc -<? cosl1 WW + sinh ^p) 2 G p (k-q,u- e£ p ) 



+ ^2" cos2 2 * Yl Y^+ q > cosh <p QiP sinh ip q >,± p ± Jk+q cosh <p q >,± p sinh cp qiP ) 2 G ± (k -q-q f ,uj - e L qp - e^± p ) 



zH 2 



q,q' ±,p 



+ ^2^2^k-qCOshe qyP smhe q ^± p ±j k -q> cosh Q q > i±p sinh q , p ) 2 G ± (k -q- q',u>- e£ p - e^ ±p ) (51) 



q,q' ±>P 



which depends on the exciton propagator and the Bo- 
golyubov angles derived in the previous section. A sim- 
ilar formula to (51) applies to E". However, it is easily 
verified that 



E~(k,uj) = E + (/c + (7t,7t),cj) 



(52) 



since 



7/c+(7T,7T) 



-7fc. In general the SCBA (51) cannot 



be solved analytically, and hence we have obtained the 
exciton spectral function 



A(k,uj) 



1 



Im [G(fc,c 



(53) 



using an iterative procedure with Monte Carlo integra- 
tion over the spin wave momenta discretized on a 32 x 32 
momentum grid. We start with E = and after approxi- 
mately 20 iterations the spectral function converged. The 
results for typical values of a, J and t are shown in Fig- 
ures E] to [131 

We start from the situation with a > 1 where the mag- 
netic background is a disorder phase with all spin singlet 
configuration in the same rung. In this case, the free dis- 
persion of the exciton with bandwidth proportional to t 
survived because all the magnetic triplet excitations are 
gapped, with an energy of Jz^/a(a — 1). For t < J, 
the exciton-magnetic interactions will barely change the 
free dispersion while for t > J such exciton-magnetic 
interactions can still occur, leading to a small 'spin po- 
laron' effect where the exciton quasiparticle (QP) peak 



T 

A(k,co) 




dispersion: bandwidth 2zt 



FIG. 8: Exciton spectral function for parameters J — t and 
a = 1.4. The only relevant feature is the strong quasiparticle 
peak with dispersion equal to 8t, where t is the hopping en- 
ergy of the exciton. The horizontal axis describes energy, the 
vertical axis is the spectral function in arbitrary units. 



is diminished and spectral weight is transferred to a po- 
laronic bump at a higher energy than the quasiparticle 
peak. For most values of t/J this effect is however neg- 
ligible already for a just above the critical point. The 
exciton spectral function for t = J and a = 1.4 can be 
seen in Figure [8] 

As a decreases towards the quantum critical point at 
a = 1, the gap of the triplet excitations also decreases. 
The effect of the exciton-magnetic interactions become 
more significant, which leads to an increasing transfer of 
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A(k,co) 



J = 0.2 t = 0.075 eV 
t = 0.375 eV 
a = 1.0 




T 

A(k,co) 




FIG. 9: Exciton spectral function at the quantum critical 
point, for J = 0.2£ and a = 1. No distinct quasiparticle 
peak is observable, and at all momenta a broad critical bump 
appears in the spectrum. 



FIG. 10: Exciton spectral function for parameters J = lOt 
and a = 0. A small nondispersive peak, indicative of lo- 
calization, dominates the spectrum. Spin wave interactions 
generate a second smaller peak at an energy Jz above the 
main peak. 



spectral weight from the free coherent peak to the inco- 
herent parts. When a hits the quantum critical point the 
gap of all spin excitations vanishes. There the motion of 
the exciton is strongly scattered by the spin excitations 
which completely destroy the coherent peak and leads to 
an incoherent critical hump in the spectrum as shown in 
Figure [9] For a further decreases to values a < 1, mag- 
netic background becomes antiferromagnetically ordered 
with two gapless transverse modes and one gapped lon- 
gitudinal mode. In this case, the motion of the exciton is 
still strongly scattered with the spin excitations leaving 
a footprints in the exciton spectrum. 

Most striking thing happens then at a = 0, when the 
two layers are effectively decoupled and we would expect 
similar behavior for an interlayer exciton as for a hole 
or electron in a single layer. Ind eed conform with the 
single hole in the t — J model 2 ^^ we find that a moving 
exciton causes spin frustration with an energy propor- 
tional to J. In the limit where J ^> t the kinetic energy 
of the exciton is too small to be able to move through 
magnetic background. Therefore, we expect a localiza- 
tion of the exciton which is reflected in spectral data by 
an almost non-dispersive quasiparticle peak. This peak 
has a bandwidth proportional tot 2 /J and carries most of 
the spectral weight, 1 — G(t 2 /J 2 ). The remaining spec- 
tral weight is carried by a second peak, at an energy Jz 
above the main peak. Numerical results are shown in 
Figure [lOj 

More complex behavior at a = arises in the anti- 
adiabatic limit t ^> J, where the kinetic energy of the 
exciton is large compared to the energy required to excite 
(and absorb) spin waves. Consequently, many spin waves 
are excited as the exciton moves and the exciton becomes 
'overdressed' with multiple spin waves. At nonzero J 
however, a very small quasiparticle peak remains with 
a bandwidth of order J. Nonetheless the majority of 
spectral weight is carried in the incoherent many-spin 



wave part, as is shown in Figure 11 



T 

A(k,co) 



J = 0.2 t = 0.05 eV 
t = 0.25 eV 
a = 



non-dispersive 
Ising peaks 




incoherent many-spin wave part, width 8zt 



FIG. 11: Exciton spectral function for parameters J — 0.2t 
and a = 0. Due to many-spin wave interactions the spectrum 
is barely dispersive and contains a large incoherent bump. 



of t/J. What happens here? A simple extrapolation 
of the two a = limits yields that the bandwidth of 
the quasiparticle peak will reach its maximum value at 
J « t. Similar extrapolations suggest that about half of 
the spectral weight will be carried by the QP peak. How- 
ever, inclusion of a finite value of a is not so trivial. Our 



Finally, realistic physical systems are expected to have 
a small nonzero value of a and an intermediate value 



quantitative computations in Figure 12 show a remark- 
able emergent Ising-like ladder spectrum, which hints at 
confinement of the exciton. Physically, this can be un- 
derstood as follows. In the a = limit, the magnetic 
interactions are dominated by the transverse excitations 
which are just single layer spin waves. The (interlayer) 
longitudinal spin waves become only relevant at finite a. 
Since these longitudinal spin waves are gapped and non- 
dispersive, interactions with longitudinal spin waves lead 
to a ladder spectrum for the exciton. In first order in a, 
the interaction with transversal spin waves is reduced and 
the longitudinal interactions become dominant. The re- 
sulting Ising-like spectrum contains multiple peaks, and 
the energy distance between the two lowest peaks scales 
as t(J/t) 2 / 3 . The spectral weight carried by higher order 
peaks vanishes as t/J —> 0. This is exactly the behavior 
seen in Figure [12] 
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A(k f Cfl) 



J = 0.2 t = 0.075 eV 
t = 0.375 eV 
a = 0.2 



almost non-dispersive 
Ising ladder spectrum 




FIG. 12: Exciton spectral function for J — 0.2t and a — 
0.2. On top of the incoherent bump a strong ladder spectrum 
has developed, signaling Ising confinement. The exact Ising 
ladder spectrum is shown in green dotted lines. The Ising 
peaks are very weakly dispersive, with bandwidth of order J. 



electric function or any other charge-excitation measure- 
ments. One particular example of the former is electron 
energy loss spectroscopy (EELS) which showed earlier 
clear signa tures of the in-plane charge transfer excitons 
in cupratei^l^J The EELS cross-section is directly re- 
lated to the dielectric function 43 via the dynamic struc- 
ture factor S(q,uj), 



da oc —S(q,uj) oc -rrlm 



-l 



where the dynamic structure factor equals 

A i 

x(A|^e^> ) 



(56) 



S(q,0j) 



e -ig-r ie i(u-H)t^ 



(57) 



The ladder-like spectrum for small nonzero values of 
a is quite remarkable, and is a typical bilayer feature. 
To understand this notice that the interactions with the 
longitudinal spin wave are responsible, in first order in 
a, for the Ising spectrum. Take the SCBA equation (51), 



neglect the transversal vertex and expand the self-energy 
up to first order in a. Only the single spin wave diagram 
contributes and it equals 



E + (fc,o;) 



N ^ 



rt_ q G*(k-q,w-Jz) (54) 



from which we deduce, observing that S~ = S + and 
shifting the momentum summation, that the self-energy 
must be momentum-independent and given by the self- 
consistent equation 



S(w) 



Jz — E(cj — Jz) 



(55) 



This self-energy is exactly the same as the self- 
energy of a single dopant moving through an Ising 
antiferromagnet! 2 ^. In fact, in any system where a mov- 
ing particle must excite a gapped and flat mode an Ising 
confinement spectrum will appear. 

Of course the exciton ladder spectrum in Figure [12] is 
not exactly sharp. By the above analysis, we can infer 
that the incoherent broadening of peaks is due to in- 
teractions with the transversal spin waves. Indeed, the 
transversal spin waves can be viewed as the equivalent 
of the single layer spin waves. Therefore for small a the 
effect of transversal spin waves is to reproduce the results 
for a single hole in the t — J model, which is quasiparticle 
peak broadening. 



V. RELATION TO EXPERIMENT 

The formation of bound exciton states can be exper- 
imentally verified in indirect measurements of the di- 



where the sum A runs over all intermediate states, |^o) 
is the . We use the dipole expansion such that 



a iqri 



1 + iq • fi + . . . 



(58) 



where the electron position operator can be expanded in 
terms of the possible electron wave functions in the tight 
binding approximation, 



^r i = ^c\ a c jcT {(t) i \r\(t) j ) 



(59) 



where \<pi) are the Wannier wave functions of the electron 
on site i. The z component of {(j>i\r\(j)j) is proportional to 
the interlayer hopping energy t_\_, which in turn is equal 
to the the creation operator of an exciton, 



oc 



(60) 
(61) 



We recognize the Fourier transform of the k = excitonic 
state, so that we find 



S(q z ,oj) oc (q 



" t±)2 J ^ e " e|t 'E^ol^=o e*-*>*|A) 



(62) 



We have introduced the term e _e ' £ l to ensure convergence 
of the integral so that we can integrate over t. We find 
that the dynamic structure factor is directly related to 
the exciton spectral function 



S(q z ,w) oc (q z t ± ) 2 (iP \E k=0 



k uj — H + ie 
(x (q z t ± YA(k = 0,uj) (63) 



uj — H — ie 

\2 
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FIG. 13: Expected exciton spectral function for the c-axis 
charge- transfer exciton in YBCO bilayers. We used model 
parameters J = 0.125 eV, t = 0.1 eV and a = 0.04. The 
exciton quasiparticle peak has a dispersion with bandwidth 
t 2 /J, and the quasiparticle peak is the most pronounced at 
the line between (ir, 0) and (0,7r). Following at a distance 
of zt(J/t) 2//3 , a secondary peak develops as a sign of Ising 
confinement. 



or in other words 



Im [e-^q^uj)] ~ (t 1 _) 2 A(k = 0,u). 



(64) 



Consequently, one expects that the bound exciton states 
to show up in EELS measurements when probing the z- 
axis excitations. In addition to the bound exciton states, 
a broad electron-hole continuum will show up at high 
energies. 

When comparing the dielectric function with the com- 
puted spectral functions in Figures |8p3} bear in mind 
that the latter are shifted over the energy Eq required to 
excite an interlayer exciton. This energy is typically of 
the order of the interlayer Coulomb energy. 

Another possible way to detect interlayer excitons is to 
use optical probes. The optical conductivit y cr (q,uj) of a 
material is related to the dielectric function 44 by 



i — V c (q)cr(q,uj), 



(65) 



where V c (q) is the Fourier transform of the Coulomb po- 
tential — | jt . The real part of the c-axis optical con- 



eo\r — r' \ ' 

ductivity is therefore proportional to the exciton spectral 
function. Similar considerations hold when one measures 
the Resonant Inelastic X-ray Scattering (RIXS) 45 spec- 
trum. 

How would then the exciton spectrum look like 
for a realistic material, such as the bilayer cuprate 
YBa2Cu307_5 (YBCO)? Following earlier neutron scat- 
tering experiments^ 4 ^ one can deduce that the effective 
J\\ = 125 ± 5 meV and J± = 11 ± 2 meV, which corre- 
sponds to an effective value of a = 0.04a c where a c is 
the critical value of a. 25 . The question remains what a 
realistic estimate of the exciton binding energy is. The 
planar excitons are known to be strongly bound 42 with 
binding energy of the order of 1-2 eV. Since the Coulomb 
repulsion scales as V ~ (er) _1 , we can relate the bind- 
ing energy of the interlayer excitons to that of the planar 



excitons. The distance between the layers is about twice 
the in-plane distance between nearest neighbor copper 
and oxygen atoms, but simultaneously we expect the di- 
electric constant e c along the c-axis to be smaller than 
e a 5 due to the anisotropy in the screening. Combining 
these two effects, we consider it a reasonable assumption 
that the interlayer exciton binding energy is comparable 
to the in-plane binding energy. The hopping energy for 
electrons is approximately t e = 0.4 eV which yields, to- 
gether with a Coulomb repulsion estimate of V ~ 1.5 eV, 
an effective exciton hopping energy of t ~ 0.1 eV. Note 
that these estimates of V jt justify our use of the strong 
coupling limit in section PH 

The spectral function corresponding to these parame- 
ters is shown in Figure [13] Since t ~ J the ladder spec- 
trum is strongly suppressed compared to the aforemen- 
tioned anti-adiabatic limit. However, the Ising confine- 
ment still shows its signature in a small 'second ladder 
peak' at 0.4 eV energy above the exciton quasiparticle 
peak. To the best of our knowledge and to our surprise, 
the c-axis optical conductivity of YBCO has not been 
measured before. Detection of this second ladder peak 
in future experiments would suggest that indeed the in- 
terlayer excitons in cuprates are frustrated by the spin 
texture. 



VI. CONCLUSION 

Using a rung linear spin wave theory for the bilayer 
Heisenberg model we constructed a theory of strongly 
bound excitons in a strongly correlated bilayer system. 
Surprisingly, for small but finite a = the exciton be- 
comes confined in a fashion similar to Ising confinement. 
The resulting ladder spectrum should be visible in mea- 
surements of the dielectric function, such as EELS, RIXS 
or optical conductivity. 

Possible candidate materials are for example het- 
erostructures of n and p-type doped cuprates such 
as Nd 2 -^Ce :c Cu04/La2-xSr :c Cu04. In YBCO or 
Bi2Sr2CaCu208+(5, the copperoxide layers come in pairs 
which suggests the possibility of interlayer charge- 
transfer excitons. A spectrum of c-axis excitons in un- 



doped YBCO is shown in Figure 13 



Our model can be extended to different stacking struc- 
tures. For example, in 214 compounds the sites in adja- 
cent cuprate layers do not lie above each other, and we 
might need to include new interlayer magnetic interac- 
tions such as the Moriya-Dzyaloshinskii interaction. Dif- 
ferent lattice structures can also be studied, of which the 
hexagonal lattice (as in graphene) is the most relevant. 

One may wonder to what extent the used approxima- 
tions are generally valid, such as the linear and ideal 
spin wave approximation. For the single layer Heisen- 
berg model, it was shown that the next-to-leading order 
corrections where indeed significantly smallei^, justify- 
ing the use of both approximations in that case. To- 
gether with the fact that we were able to reproduce the 
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known phase diagram and excitation spectrum, this sug- 
gests our approach for the bilayer Heisenberg model is 
justifiable. Nevertheless, an exact computation of the 
next-to-leading order corrections can quantify the errors 
of the used spin wave approximations. 

Another approximation we used was the expansion in 
large V, the exciton coupling strength. This coupling 
originates in the interlayer Coulomb interaction, from 
which we only consider the on-site and nearest neighbor 
terms. Therefore our model cannot describe accurately 
the process of how excitons are formed out of separate 
doublons and holons. We think this is a very interesting 
open question, especially at finite temperatures. In ad- 
dition, the formation process is also accompanied by an 
exciton annihilation process which we neglected in our 
current work. 

Besides the interesting properties of the exciton forma- 
tion process, we think that further research should be di- 
rected towards finite densities of excitons^. The dynam- 
ical spin-hole frustration effects that are well known in 
the context of doped Mott insulators occur in a strongly 
amplified form dealing with interlayer excitons in Mott- 
insulating bilayer systems. This gives further impetus to 
the pursuit to create such finite density correlated exci- 
ton systems in the laboratory. One can wonder whether 
such physics is already at work in the four-layer material 
Ba2CasCu408F2 where self-doping effects occur creating 
simultaneously p and n-doped layers^. Much effort has 
been devoted to create equilibrium finite exciton densities 
using conventional semiconductors^, while exciton con- 
densation has been demonstrated in coupled semiconduc- 
tor 2DEGs^. In strongly correlated heterostructures, 
however, formation of finite exciton densities is still far 
from achieved, although recent developments on oxide in- 
terfaces indicate exciting potential (see for example^). 
Besides the closely coupled p- and n-doped conduct- 
ing interface- layers in these SrTi03-LaA103-SrTi03 het- 
erostructures, further candidates would be closely cou- 
pled p- and n-doped cuprates, such as YE^CusOr-^ 
or La2- x Sr :c Cu04 with iNk^-^Ce^CuCU. The feasibility 
of this has already been experimentally demonstrated, 
e.g. in 50 , but the exact interface effects need to be in- 
vestigated i n mor e detail, both experimentally as well as 
theoretical!}^ 11 ^. 



Appendix A: Large S limit bilayer Heisenberg model 



In this appendix we will prove equation (26). The 



mean field Hamiltonian (25) depends on the antiferro- 
magnetic (AF) order parameter fh. We must find the 



ground state energy of ( 25 ) as a function of m and then 



minimize with respect to m, thus yielding the mean field 
value of the AF order parameter. 

However, since we are only interested in the critical 
value a c where fh changes from nonzero to zero, we can 
proceed as follows. In the singlet phase (fh = 0) the 
mean field Hamiltonian is reduced to 



(0) = J ± Sx • S 2 



(Al) 



which has as ground state the singlet |0 0) and as first 
excited state the triplet 1 1 0) with energy difference E\ — 
Eo — J±- We will treat the Hamiltonian terms that 
depend on fh as a perturbation, and compute the ground 
state energy in second order perturbation theory for small 
fh. If the ground state energy decreases with nonzero m, 
then there is an instability towards antiferromagnetism. 
The perturbation Hamiltonian is 



1 



Jzm 



1 



Jzm{-iyS z 



(A2) 



and the first and second order corrections to the ground 
state energy are 



Now contains one term that is just an identity oper- 
ator, and the S z operator can only change the total spin 
number s by one single unit. This means that the former 
expression yields 



(i) 



(2) 



-Jzm 1 
4 

Jzm 2 
4a 



(Jzffh 2 



4J± 

a- |(10|^|0 0)| 2 



10|^|0 0)| 2 

(A4) 



We see that whenever a > |(1 0|*S^|0 0)| 2 , the ground 
state energy always increases when fh is nonzero. Hence 
the critical value of a is given by 



a c = |(10|^|0 0)| 2 . 



(A5) 



The right hand side can be evaluated explicitly using 
Clebsch-Gordan coefficients, since 



(1 0|S*|0 0) 



E 



9m /"'SS'l (~tSS0 



2S 
2 

71 



+ 1 y 5(5 + 1) 
\/5(5 + l) 



E 



m 



from which we indeed conclude that 



a c = -5(5+1). 



(A6) 



(A7) 



Appendix B: Bilayer Heisenberg Hamiltonian in 
terms of e, b operators 

The bilayer Heisenberg operators (total spin and spin 
difference) can be expressed in terms of the local spin 
excitations and b\ by 
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St = 



qz _ 



b\ ia b +ia - bl ia b- ia 

- sin x(b+ ia + b-ia) + cos x{b\ ia e ia + el a b- icr )j 

(-1)"* ^sin2 X (l -Y, b ±iab±ia - 2eJ CT e i<T ) - cos2 X (eJ (7 + 

V^i-iy- (cos X (^ i<T - b-ia) + sm X (bl ia e ia - e\ a b_ ia ) 



(Bl) 
(B2) 

(B3) 

(B4) 



where a represents the sign of the sublattice of site i. Consequently, the bilayer Heisenberg model in terms of these 
new operators reads (with a = and a = A, B denotes the sublattice index), 

H = -JzN(—a — 2Xa cos 2x — sin 2 2%) + ^ Jz y^(cos 2\ — &X) sin 2x(e]~ cr + ej a ) 

i 

+Jz^(sin 2 2x + a\cos2x)el a e ia - - J ^ cos 2 2x(e\ A + e iA ){e\ + ^ B + e i+5 , B ) 

i ieA,5 

+ 7) Jz ^2( a + sin2 2 X + cos2x)^La^±^ 

+ h.c) 



ieA,S 



+C(6 t 6 t e + e t 66) + O ([e f e + 6 f 6] 2 ) 



(B5) 



We explicitly neglect the interaction terms, which are cu- layer Heisenberg modePS. However, in the bilayer model 
bic and quartic in the spin wave operators. The above there will be a competition between the ordered phase 



Hamiltonian contains a constant term (depends only on 
a, x an d A) that describes the ground state energy com- 
petition between the singlet and antiferromagnetic phase, 
see appendix[C] The term linear in spin operators gives us 
the self-consistent condition for x- The quadratic terms 
will be diagonalized using the Fourier and Bogolyubov 
transformation as described in the main text. 



(cos2x = otX) and the disordered phase (x = 0). Note 
for aX > 1 we automatically end up in the disordered 
phase. 



Appendix C: Quantum phase transition 



The ideal spin wave approximation introduces a shift 
in the ground state energy, similar to that in the single 



For aX < 1, the ground state energy of both phases is 
given by the expression 



E = - J 'zN (-a -2aX cos 2% - sin 2 2%) 



4 

k 



1 



(sin 2x + 06 X cos 2%) sinh y?k — — cos 2x7fc(cosh 2cpk — sinh 2(fk) 



(a + sin 2 2x + oiX cos 2%) sinh 2 Ok — - Xjk (cos 2x cosh 20 k + sinh 20k ) 



(Cl) 



where we have to fill in the right values of x> Ok an d (fk 



depending on the phase. As can be seen in figure [l4j 
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the spin waves drive the system earlier into the singlet 
phase, namely at a c « 0.605. For smaller values of A this 
critical value increases, proportional to A -1 . 

The critical value a c = 0.605 for our spin wave the- 
ory closely resembles the numerical results of a c = 0.63. 
Since this ground state energy competition the system 
is driven into the disordered state for a different a than 
mean field theory suggests, we should replace bare val- 
ues of a — by the renormalized a* = a/a c when 
computing the exciton spectral function. 



FIG. 14: Ground state energies of the bilayer Heisenberg 
model following equation ( CI ). Shown is the energy of the an- 



tiferromagnetic phase (in red) and the singlet phase (in green) 
for the isotropic A = 1 model. The energies are measured in 
units of JzN. At a « 0.605 there is a phase transition from 
the AF to the singlet phase. 



Appendix D: Explicit expressions for exciton-spin 
wave interactions 



We can rewrite the Hopping Hamiltonian (16) from 
the singlet-triplet basis using the local spin excitation 
operators defined in equations (28)- (30), 



H t =tY,E)Ei I cos2 X (l - e\e 3 ) + sin2 X (et + e 3 ) - J2 h l h j* ) ( D1 ) 

(ij) V / 

where a is the sum over spins ±1 and < ij > denotes nearest neighbor pairs. We need to rewrite this in terms of the 
longitudinal (£) and transversal (a and fi) modes derived in the main text. Therefore we first split all operators into 
the ones that live on sublattice A and the ones that live on B, 



H t = t^2 E i+5,B E iA I cos2x(l - et A e i+5)jB ) + sin2x(et A + e i+8iB ) - XX 

ieA,5 V cr 



A a b i+5 ,B,a + h.c. (D2) 



As described in the main text, the Fourier transform for the sublattice operators is E^ A = y j^^2 k E\ A e lkri and we 

introduce the in-phase p = 1 and out-phase p = — 1 exciton operators p — -^{E\ A + pE^ B )\ similar expressions 
hold for the spinon operators. The hopping Hamiltonian now can be written as 

#o = ztcos2 X J2P^El p E k:P (D3) 

k,p 

zt 



k,q p,p' 



zt 

H 2 = --^cos2 X J2J2 p P'^ k - k '+i E k+ q ,P E k,±P e l- q ,±p' e k', P ' (D5) 

k,k',qp,p' ± 

z~b 

H * = -ff Yl J2J2J2PP'7 k - k '+* E k+ q ,p E k,±P b l-q,±p>,* b ki,pi,a (D6) 
k,k f ,q cr p^p' ± 

Note that this Hamiltonian contains four different type of processes. The first line Hq contains a free part of the 
exciton motion. The bandwidth of the free exciton dispersion increases linearly in a in the antiferromagnetic phase 
until it saturates at 2zt in the disordered phase. The next term Hi describes the creation and annihilation of a single 
longitudinal mode due to exciton motion. This term is only present in the antiferromagnetic phase and is comparable 
to the hole-spin vertex in the single layer t — J model. Finally, there are two H2 interactions where an exciton 
scatters off a transversal (H2) or longitudinal (H2) mode. These processes can also be changed into the creation or 
annihilation of a pair of spin modes. All processes can be characterized by a conservation of total phase index p and 
conservation of total momentum. 
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The remaining step is to write out the interaction vertices explicitly in terms of the Bogolyubov transformed spin 
waves. The single-magnon process equals 



zt sin 2\ 



N 



ki...k 3 P1-.-P3 i=l 

><Pi(P3lk 2 cosh cp k3P3 + j kl sinh<p fe3P3 ) + h.c. 
The process that involves two longitudinal spin waves is given by 



(D7) 



Hi 



zt cos 2x 
N 



E E si2) ^ - k * + ** - wiift - !) 4 1 p 1 ^ 2 p 2 ci 3 p 3 ^ 

fcl .../C4 Pi ...P4 2 = 1 



xpi (P47fei-fe 4 cosh (^^3 cosh(^ fe4P4 + p 3 7fci-fc 3 smn( /^ 3 p 3 sinh(/? fe4P4 ) 
cos 2% 1 



TV 



E E 5(2) ^ - k * + ** + wllft - !) ^L^^cLsCl 

fcl ...fe4 Pi ..-P4 * = 1 

xpi (P47fe!+/c 4 cosh(/p fe3P 3 sinh(/p fe4P4 + £>37fci+/c 3 sinhy? fe3 P3 coshy? fe4 p 4 ) + ^- c -- 



4P4 



(D8) 



Finally, we can also write out the Hamiltonian for the interaction vertex with the transversal spin waves. We can 
write this Hamiltonian term explicitly using phase and momentum conservation, 



-f E E * (2) ( fc i " fc 2 + h - h)5(fl Pl - 1) El ip E k2P2 (al p a k4P4 + $l P3 fa iPi ) 

ki...k 4 pi-.-PA i=l 

xpi (P47fe!-/c 4 cosh (9^3 cosh6> fe4P4 + p 3 7fci-/c 3 sinh/9 fe4P4 sinhfl^) 

4 

£ £ J (2) (A . 1 _ A!2 + A . 3 + W JJp._ 1) 4 ipi ^ 2P2 at 3p3 ^ 4p4 

fcl...fe4 P1-P4 1=1 

xpi (P47fe!+/c 4 cosh 6^3 sinh0fe 4P4 + P37/ci+fe 3 cosh# fe4P4 sinh0 fe3 p 3 ) + ^- c -- 



zt 
"TV 



(D9) 



r 



These expressions are used to transform the Feynman Hilgenkamp. The authors wish to thank Jeroen van der 
diagrammatic representation of the SCBA of Figure [7] Brink, Sergei Mukhin and Matthias Vojta for helpful dis- 



into the explicit formula (51) 



cussions. 
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